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I summarize the discovery potential for physics beyond the electroweak scale at the 
Pierre Auger Observatory. This observatory is designed to study ultra-high energy 
cosmic rays with unprecedented precision, with the primary goal of shedding light 
on their composition and origins. In addition, since the center-of-mass energies of 
Auger events are well beyond those reached at terrestrial colliders, they provide 
an opportunity to search for new physics. I discuss here some of the relevant 
observables and techniques which may bo used to weed out theories beyond the 
standard model. 



1 GZK end of the cosmic ray spectrum? 

Shortly after the cosmic microwave background (CMB) was discovered, 
Greisen, Zatsepin, and Kuzmin (GZK) pointed out that the rehc photons 
make the universe opaque to cosmic rays (CRs) of sufficiently high energyP 
This occurs, for example, for protons with energies beyond the photopion 
production threshold, 
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where rrip (to^) denotes the proton (pion) mass and Ecmb ^ 10^'^ eV is a 
typical CMB photon energy. After pion production, the proton (or perhaps, 
instead, a neutron) emerges with at least 50% of the incoming energy. This 
implies that the nucleon energy changes by an e-folding after a propagation 
distance < {<Jpjn~^y)~^ ^ 15 Mpc. Here, n-^ « 410 cm~'^ is the number 
density of the CMB photons, ap^ > 0.1 mb is the photopion production cross 
section, and y is the average energy fraction (in the laboratory system) lost 
by a nucleon per interaction. For heavy nuclei, the giant dipole resonance can 
be excited at similar total energies and hence, for example, iron nuclei do not 
survive fragmentation over comparable distances. Additionally, the survival 
probability for extremely high energy (« 10^^ GeV) 7-rays (propagating on 
magnetic fields » 10~^^ G) to a distance d, p{> d) « exp[— d/6.6 Mpc], be- 
comes less than 10^'* after traversing a distance of 50 Mpc. All in all, as otir 
horizon shrinks dramatically for energies > 10^^ GeV, one would expect a 
sudden cutoff in the energy spectrum if the CR sources follow a cosmological 
distribution.El 

At the beginning of summer 2002, in a pioneering paper Bahcall and 
WaxmarP noted that the energy spectra of CRs reported by the AGASA, the 
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Fly's Eye, the Haverah Park, the HiRes, and the Yakutsk Collaborations are 
consistent with the expected GZK cutoff. As one can readily see in Fig.^ after 
small adjustments (within the known uncertainties) of the absolute energy 
scale, all these spectra are shown to be in agreement with each other for 
energies below 10^^ GeV. A point worth noting at this juncture: the analysis 
that follows, which is based on counting events above roughly the expected 
cutoff, takes the data at face value and consequently does not attempt to 
evaluate possible correlated systematic errors. The particulars of the present 
sensitivity to a super-GZK flux (i.e., CR intensity beyond 10^^ GeV) are given 
in Table 1 . Armed with the expected number of events and assuming Poisson 
statistics, it is easily seen that the existing data show significant evidence for 
a supression in the CR flux beyond 10^^ GeV. Such a supression is found 
to be a ~ 3.5a effect according to the Fly's Eye normalization, increasing 
up to ~ 8(7 if the selected normalization is that of Yakutsk. It is important 
to emphasise that with data sets as small as these and with the inherent 
uncertainty associated with energy fluctuations, this result should be treated 
with great caution. 

Systematic uncertainties in normalization have been thus far a headache 
common to all surface arrays. Already in 1986, the SUGAR Collaboration 
instead of giving a unique primary energy for the observed events, they re- 
ported the showers' equivalent vertical muon number A^^v together with two 
possible conversions from A^^v to primary energyPl Interestingly, if one adopts 
the Hillas conversion factor, for which the SUGAR integral e nerg y spectrum 
is in good agreement with the AGASA spectrum at 10^ GeV,'^the number 
of observed super-GZK events is consistent with the one expected from an 
extrapolation of the sub-GZK spectrum (cx i?^^ '') at the la level. On the 
other hand, if one adopts the Sidney normalization there are no events with 
energy > 10" GeV. It should be noted that the SUGAR exposure given m 
Table 1 was obtained on the basis that the super-GZK detection probability 
over the entire array (with maximum collecting area of 70 km2) is 85%.li I 
have also assumed that the experiment operated in stable mode from January 
1968 until February 1979, yielding a total area-time product of about 775 km^ 
yr. However, since the sensitive area varied as the array was developed and as 
detectors require maintenance, the expected numbers of events given in the 
last row of Table 1 should be taken as upper limits. 

To make the situation more confusing, the cali brations from the various 
experiments are themselves moving targets. 1^131^3 Even in the time between 
this conference and the deadline for its proceedings the Yakutsk spectrum has 
change substantially. Nowadays it seems easier to predict the stock market 
or even the weather in Boston than the next revision in the energy spectrum. 
There is no doubt that more and better data is needed. 

The arrival directions of the super-GZK events have no apparent counter- 
parts in the Galactic plane nor in the Local Supercluster. Furthermore, the 
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Figure 1. The left panel shows the upper end of the energy spectrum as reported by the 
different collaborations. The right panel shows the data after adjusting the absolute energy 
calibrations of the various experiments so as to bring the results from the different experi- 
ments into agreement at 10^" GeV. Here, the Fly's Eye energy is adopted as the standard. 
The fractional shifts in absolute energy scale, are well within the published systematic 
errors 



data is consistent with an isotropic distribution of sources, in sharp contrast 
to the anisotropic distribution of hght within 50 Mpc. The apparent isotropy 
in the arrival directions can be explained if the particle orbits are bent in ex- 
tragalactic magnetic fields .El However, the low statistics leaves a window open 
for less mundane explanations, which I discuss next. 



Table 1. Numbers of events observed with average energy > 10^^ GeV and incident zenith 
angle < 45°. The super-GZK exposure of Volcano Ranclj^ Haverah ParJ^ Fly's 'E.yW^ 
Yakutsdnj AGAS^ HiReJH and SUGAR (see main text) is given in the J,^'^ column. 
The last column indicates the expected number of events calculated based on the spectrum 
observed above 10^" GeV from Fly's Eye and Yakutsk, and assuming no change in the 
spectral index. 



Experiment 


Events observed 


Exposure [m*^ s sr] 


Expected events 

Fly'sBye Yakutsk 


Volcano Ranch 


1 


2.0 X 10^=^ 


0.4 


1.0 


Haverah Park 





5.6 X 10^5 


1.2 


2.9 


Fly's Eye 


1 


2.6 X 10^'^ 


5.4 


13.4 


Yakutsk 


1 


2.8 X IQi'^' 


5.8 


14.4 


AGASA 


11 


5.1 X 10^'' 


10.6 


26.2 


HiRcs 


2 


6.9 X 10^^ 


14.3 


35.5 


SUGAR 


5 


< 3.1 X W'^ 


< 6.4 


< 15.9 
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2 The usual suspects 



Physics from the most favored theories beyond the standard model (SM) hke 
string/M theory, supersymmetry (SUSY), grand unified theories (GUTs), and 
TeV-scale gravity have been invoked to explain the super-GZK events. The 
conjectured origins fall into two basic categories: (i) exotic sources clustered 
nearby, (ii) GZK-evading messengers. 

The top suspect belongs to the first category. In this scenario charged 
and neutral primaries simply arise in the quantum mechanical decay of su- 
permassive elementary X particles. Sources of these exotic particles could be: 

(i) topological defects (TDs) left over from early universe phase transitions 
associated with the spontaneous symmetry breaking that underlies GUTs^^ 

(ii) some long-lived metastable super-heavy relic (MSR) particles which may 
constitute (a fraction of) the dark matter in galactic haloes. Arguably, the 
observed magnitude of the CMB fluctuations fixes the reheat temperature 
following inflation to 10^^^^ GeV, allowing gravitational and thermal produc- 
tion of MSRs during the inflationary stage of the universe with just the right 
mass (m x > 10^^ GeV) for producing 10"'^^ GeV secondary particles via 
decay.l^ Due to their topological stability, the TDs (magnetic monopoles, 
cosmic strings, domain walls, etc.) can survive forever with X particles 
{mx ~ 10^^ — 10^^ GeV) trapped inside them. Nevertheless, from time to 
time, TDs can be destroyed through collapse, annihilation, or other processes, 
and the energy stored would be released in the form of massive quanta that 
would typically decay into quarks and leptons. Discrete gauge symmetries or 
hidden sectors are introduced to stabilize the MSRs and so generally higher 
dimensional operators are required to break the new symmetry super-softly 
to maintain an appreciable decay rate today (coUisional annihilation has been 
considered, too). MSRs would also have quarks and leptons as the ultimate 
decay products. The strongly interacting quarks fragment into jets of hadrons 
containing mainly pions together with a 3% admixture of nucleons. The in- 
jection spectrum is therefore a rather hard fragmentation-type shape with an 
upper limit usually fixed by the GUT scale. Of course, the precise decay 
modes of the X particles and the detailed dynamics of the first second ary 
particles depend on the exact nature of the particles under consideration.!^ 
However, one expects the bulk flow of outgoing particles to be almost inde- 
pendent of such details: all top-down scenarios predict a spectrum dominated 
by photons and neutrinos produced via pion decay. 

The neutrino is the only known stable particle immune to the GZK degra- 
dation. The corresponding vi/ annihilation mean free path on the cosmic neu- 
trino background, = {ui, a^p)^^ « 4 x 10^^ cm, is just above the present 
size of the horizon (recall that H^^ ~ 10^^ cm). One may then entertain the 
notion that neutrinos are indeed the super-GZK primaries. Unfortunately, 
perhaps, a^N is, within the SM, about five orders of magnitude too small 
to explain the observed atmospheric cascades. On the other hand, the limit 
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imposed by unitarity is relatively weak and consequently does not impact on 
new interactions beyond the elec troweak scale to increase significantly the 
neutrino-nuclcon cross section.EI 

On a different track, if some flavor of neutrinos has a mass ~ 0.1 eV, the 
relic neutrino background is a target for everyday weakly interacting neutri- 
nos to form a Z-boson that su bseq uently decays producing a "local" flux of 
nucleons, photons and ncutrinosM^ To reproduce the observed spectrum, the 
Z-burst mechanism requires very luminous sources of extremely high energy 
neutrinos throughout the universe (see Fig. [SJ.^^ 

A novel beyond-SM-model proposal to break the GZK barrier is to as- 
sume that ultrahigh energy CRs are not known part icles but a new species of 
particle, generally referred to as the uhecron, The meager information 

we have about super-GZK particles allows a naive description of the prop- 
erties of the U. The muonic content in the atmospheric cascades suggests 
U's should interact strongly. At the same time, if U's are produced at cos- 
mological distances, they must be stable, or at least remarkably long lived, 
with mean-lifetime r > 10^ (mjy/S GeV) (c?/Gpc)s, where d is the distance 
to the source and mu, the uhecron's mass. Additionally, since the thresh- 
old energy increases linearly with mu, to avoid photopion production on the 
CMB mu > 1.5 GeV. In this direction, light supersymmetric baryons (made 
from a light gluino -I- the usual quarks and gluons, mij < 3 GeV) produce 
atmospheric cascades very similar to those initiated by protons .1^2 

Another interesting possibility in which super-GZK CRs can reach us from 
very distance sources may arise out of photons that mix with light axions in 
extragalactic magnetic flelds.'^Sl These axions would be sufficiently weakly 
coupled to travel large distances unhindered through space, and so they can 
convert back into high energy photons close to the Earth. An even more 
radical proposal postulates a tiny violation of local Lorentz invariance, such 
that some processes become kinematically forbidden.l^In particular, photon- 
photon pair production and photopion production may be affected by Lorentz 
invariance violation. Hence, the absence of the GZK-cutoff would result from 
the fact that the threshold for photopion production "disappears" and the 
process becomes kinematically not allowed. 



3 Fingerprints in the sky: The distribution of arrival directions 

The distribution of arrival directions is perhaps the most helpful observable 
in yielding clues about cosmic ray origin. Neutral GZK-evading messengers 
should point back to distant active galaxies, thereby enabling point-source 
astronomy. The earliest super-GZK events did in fact point towards high- 
redshift compact, radio-loud quasars.'^Si However, the current world data set 



"^Similarly, gravi-burst fragmentation jets can contribute to the super-GZK spectruml 
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show that such association is controversialP^ Another reveahng signature of 
discrete sources would be the clustering on a small angular scale. The data 
recorded by AGASA suggest that the pairingof events on the celestial sphere 
is occurring at a higher than random rateM^ Moreover, event directions in 
a combined data sample from AGASA, Haverah Park, Yakutsk and Volcano 
Ranch also support no chance-association. " However, it is interesting to re- 
mark that to calculate a meaningful statistical significance of CR clustering, 
one must define the search procedure a priori in order to ensure it is not 
(inadvertently) devised especially to suite the particular data set after having 
studied it. In the above mentioned analyses, for instance, the angular bin size 
was not defined ahead of time. 

Surface arrays in stable operation have nearly continuous observation over 
the entire year, yielding a uniform exposure in right ascension, a. A tradi- 
tional technique to search for large-scale anisotropics is then to fit the right 
ascension distribution of events to a sine wave with period 27r/ m (m *'' har- 
monic) to determine the components (x, y) of the Rayleigh vectoil^ 

2 w ^ N 

i=l 1=1 

The m}^ harmonic amplitude of N measurements ai is given by the Rayleigh 
vector length TZ — {x^ + y^Y^'^ . The expected length of such a vector for 

values randomly sampled from a uniform phase distribution is 7^0 = 
The chance probability of obtaining an amplitude with length larger than that 
measured is p{> TZ) — e^*^" , where — TZ'^/TZl. To give a specific example, 
a vector of length ko > 6.6 would be required to claim an observation whose 
probability of arising from random fluctuation was 0.0013 (a "3cr" result). 
For the ultra-high energy (> 10^"'^ GeV) regime, all experiments'* to date 
have reported ko 6.6.1^ This does not imply an isotropic distribution, but 
it merely means that available data are too sparse to claim a statistically 
significant measurement of anisotropy by any of these experiments. In other 
words, there may exist anisotropies at a level too low to discern given existing 
statistics. For example, a clean signature of the MSR-X hypothesis is the 
anisotropy imposed by the asymmetric position of the sun in the Galactic 
halo. As seen in the Northern hemisphere, the amplitude ~ 0.3 predicted 
by isothermal haloes with realistic core radii is in agreement with existing 
data (the most restrictive being the AGASA sample with TZra=i ~ 0.3 and 
fco w 1.0)P 

The a harmonic analyses are completely blind to intensity variations 
which depend only on declination, 5. Furthermore, combining anisotropy 



''Since Fly's Eye has had a nonuniform exposure in sidereal time, 72. was computed using 
weighted showers l^iU A shower's weight depends on the hour of its sideral arrival time. The 
24 different weights are such that every time bin has the same weighted number of showers. 
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searches in a over a range of declinations could dilute the results, since sig- 
nificant but out of phase Rayleigh vectors from different declination bands 
can cancel each other out. An unambiguous interpretation of anisotropy data 
requires two ingredients: exposure to the full celestial sphere and analysis in 
terms of both celestial coordinates. Though the statistics are very limited at 
present, the analysis of data from the SUGAR and the AGASA experiments, 
taken during a 10 yr period with nearly uniform exposure to t he e ntire sky, 
shows no significant deviation from isotropy beyond lO^"'^ GeV.ES 



4 Smoking guns: Photon and neutrino fluxes 

Another telltale discriminator to be thought is the primary composition of 
the super-GZK events. Every model wherein the primaries arise from QCD 
jets (these models include Z-bursts, and decaying TDs and MSRs) produce 
many more mesons than baryons, and consequently the injection spectrum 
would be dominated by 7-rays produced through tt" decay. Additionally, any 
technique used to distinguish photon-initiated showers from hadron-initiated 
showers would be well suited to test the validity of the photon-axion mixing 
model. 

At large zenith angles, hadrons and 7-rays develop cascades in the upper 
layers of the atmosphere. The electromagnetic component, with mean inter- 
action length 45 — 60 g/cm^, is absorbed long before reaching the ground 
by the greatly enhanced atmospheric slant depth (w 3000 g/cm^ at 70° from 
the zenith). Then surface arrays are practically only sensitive to the high en- 
ergy muons created in the first few generations of particles. The shape of the 
shower front in this type of cascade is extremely flat (with radius of curvature 
above 100 km), and the particle time spread is very narrow {At < 50 ns). 
Since showers initiated by 7-rays produce fewer muons than those initiated 
by hadrons, one expects the rate of 7-ray showers detected by surface ar- 
rays to be reduced relative to the rate from hadron showers. Therefore, the 
determination of the CR-fiux through vertical shower measurements using 
fluorescence eyes (which are fairly independent of the primary composition) 
provides a powerful tool for discriminating between hadron and 7-ray showers 
when comparing with the inclined shower rate. For example, a comparison of 
the showers recorded by the Havcrah Park experiment (in the angular range 
60° < 9 < 80°) with predictions from the observed Fly's Eye spectrum yields 
strong bounds on the 7-ray flux: above 10^° GeV, less than 48% of the pri- 
mary cosmic rays can be 7-rays and above 10^*^'^ GeV less than 50% can be 
7-rays. Both of these statements are made at the 95%CL.I^ 

Even though 7-rays dominate at production, there are some viable TD 
scenarios which predict proton fluxes that are comparable or even higher than 
the 7-ray flux at all energies. Conversely, due to the lack of absorption, 
relics clustered in the Galactic halo predict compositions directly given by the 
fragmentation function. Note that mechanisms which successfully deplete the 
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7-rays (such as efficient absorption on the universal radio background) require 
an increase in the neutrino flux to maintain the overaU normahzation of the 
observed spectrum. Within the SM, the mean free path of neutrinos is larger 
than even the horizontal atmospheric depth. Neutrinos therefore interact with 
roughly equal probability at any point in the atmosphere and may initiate 
quasi-horizontal showers in the volume of air immediately above the detector. 
These will appear as hadronic vertical showers, with large electromagnetic 
components, curved fronts (a radius of curvature of a few km), and signals 
well spread over time (on the order of microseconds). 

The event rate for quasi-horizontal deep showers from ultra-high energy 
neutrinos is 

iV = ^ / TVa ^ a^N^x {Ei) £{E,) , (3) 

where the sum is over all neutrino species i = , i/g , i/^ , i/^t , z/^- , Pt- , and all 
final states X. Na — 6.022 x 10^'^ is Avogadro's number, and d^i/dEi is 
the source flux of neutrino species i, a as usual denotes the cross section, 
and £ is the exposure measured in cm'^ w.e. sr time. The Fly's Eye and 
the AGASA Collaborations have searched for quasi-horizontal showers that 
are deeply-penetrating, with depth at shower maximum ATmax > 2500 g/cm^. 
There is only 1 event that unambiguously passes this cut with 1.72 events 
expected from hadronic background, implying an upper bound of 3.5 events 
at 95%CL from neutrino fluxes.lSSlNote that if the number of events integrated 
over energy is bounded by 3.5, then it is certainly true bin by bin in energy. 
Thus, using Eq. ^ one obtains 



E 



/ dE,NA—^<J^N^x{E^)£{E,) <3.5, (4) 

J A "-t/i 



at 95% CL for some interval A. Here, the sum over X takes into account 
charge and neutral current processes. In a logarithmic interval A where a 
single power law approximation 

^a,^^xiE.)£{E,)^E^ (5) 



is valid, a straightforward calculation shows that 



I ^^^^ ^ = i^'^^^^ ^ d<S^,/dE.) ^ A , (6) 



where S — (a-|-l)A/2 and (A) denotes the quantity A evaluated at the center 
of the logarithmic interval. The parameter a — 0.363 + — 7, where 0.363 is 
the power law index of the SM neutrino cross sections, " and /3 and —7 are 
the power law indices (in the interval A) of the exposure and flux d^t/dEi, 
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respectively. Since smhS/S > 1, a conservative bound may be obtained from 
Eqs. gl) and ©: 

Na $](a,Ar^x(£^.)) {£{E..)) {E,d^,/dE,) < 3.5/A . (7) 

By taking A = 1 as a likely interval in which the single power law behavior is 
valid (this corresponds to one e-folding of energy), and setting (Eid^i/dEi) = 
^{Eydf^i,/ dE^,) {'^i, = total neutrino flux) from Eq. ((TJ it is straightforward 

to obtain 95%CL upper limits on the neutrino fluxf^ Similar concepts are 
used by the Goldstone Lunar Ultra-high energy neutrino Experiment (GL UE) 
to set upper bounds on the neutrino flux at the ultimate energy frontier!"^ 
All these bounds are collected in Fig. |21 

Unfortunately, due to the limited size of the current ultra-high energy CR 
sample, no conclusive statement can yet be made on the hypothetized models 
discussed in Sec. 2. The jury awaits further evidence. 



5 Heavenly black holes: Probes of TeV-scale gravity 

It is intriguing - and at the same time suggestive - that the observed flux of 
CRs beyond the GZK-energy is well matched by the flux predicted for cosmo- 
genic neutrinos. Of course, this is not a simple coincidence: any proton flux 
beyond E}^^^_^ is degraded in energy by photoproducing 7r° and tt^ , with the 
latter in turn decaying to produce cosmogenic neutrinos. The number of neu- 
trinos produced in the GZK chain reaction compensates for their lesser energy, 
with the result that the cosmogenic flux matches well the observed CR flux be- 
yond 10^^ GeV. Recently, the prospect of an enhanced neutrino cross section 
has been explored in the context of theories with large compact dimensions. 
In these theories, the extra spatial dimensions are responsible for the extraor- 
dinary weakness of the gravitational force, or, in other words, the extreme size 
of the Planck mass. For example, if spacetime is taken as a direct product 
of a non-compact 4-dimensional manifold and a flat spatial n-torus T" (of 
common linear size 27rrc), one obtains a definite representation of this picture 
in which the effective 4-dimensional Planck scale, Mpi ~ 10^^ GeV, is related 
to the fundamental scale of gravity, -M^d, according to Mpj = Stt M^"^" r".!^ 
Within this framework, virtual graviton exchange would disturb high energy 
neutrino interactions, and in principle, could increase the neutrino interaction 
cross section in the atmosphere by orders of magnitude beyond the SM value; 
namely u^n ^ [E^/{W°) GeV] mbP However, it is important to stress 
that a cross section of ~ 100 mb would be necessary to approach obtain- 
ing consistency with observed showers which start within the first 50 g/cm^ 
of the atmosphere. This is because Kaluza-Klein modes couple to neutral 
currents and t he s cattered neutrino carries away 90% of the incident energy 
per interaction.!^ Moreover, models which postulate strong neutrino interac- 
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Figure 2. Current limits on neutrino fluxes. Upper bound 95%CL on d<I>, 



derived by replacing the combined exposures for deeply-developing showers of AGASA 
and Fly's Eye into Eq. Q.'^^ Upper bound 90% CL on d$„ ^„ /d-E^- from the non- 
observation of electromagnetic showers at GLUE induced by neutrinos interacting in the 
moon's rim!33 The point with error bars indicates the total neutrino flux required by the 
Z-burst mechanisni21iland the dotted line indicates a typical ci^^^+c^ /'^^•^i from top down 
cascades with = 10^'^ GevC^In these models, d^^^^-^/dE^. d^^ +i>^,/dE^.. The 
diamonds indicate the sensitivity of the Pierre Auger Observatory, i.e., any flux lying above 
thes e po ints for at least one decade of energy will give more than 1 observed event per 
yearlSil 



tions at super-GZK energies also predict that moderately penetrating showers 
should be produced at lower energies, where the neutrino- nucleon cross section 
reaches a sub-hadronic size. Within TeV scale gravity cr^jv is likely to be sub- 
hadronic near the energy at which the cosmogenic neutrino flux peaks, and so 
moderately penetrating showers should be copiously produced. Certainly, the 
absence of moderately penetrating showers in the CR data sample should be 
understood as a seriou s ob jection to the hypothesis of neutrino progenitors of 
the super-GZK events.^^ 

Large extra dimensions still may lead to significant increases in the neu- 
trino cross section. Should we be so lucky that this scenario is true, we 
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Figure 3. Left: Existing limi ts o n the fundamental Planck scale A/^i from tests of Newto n's 
law on sub-millimctcr scales(^l bounds on superno va c ooling and neutron star heatingp^l 
dielectron and diphoton pr oduc tion at the TevatronpHland non-observation of BH produc- 
tion at AGASA-I-Fly's EyeESI Right: Comparison of current bounds on M/j with future 
limits from the Southern Auger ground array, assuming 5 years of data and no excess above 
the SM neutrino background.'^'' The range in Tevatron bounds corresponds to the range of 
the brane softening parameter. The range in cosmic ray bounds is due to variations in the 
minimum BH mass. 



might hope to observe black hole (BH) production (somewhat more ma- 
sive than Mo) in e lementary particle collisions with center-of-mass ener- 
gies ^/s > TeV.ESl In particular, BHs occurring very deep in the atmo- 
sphere (revealed as intermediate states of ultrahigh energy neutrino inter- 
actions) could trigger quasi-horizontal showers and be detected by cosmic ray 
observatories. E2I Interestingly, a^N^BH oc 7Vf|j~^^^"-'^'^^"^ Therefore, using 
Eq. the non-observation of the almost guaranteed flux of cosmogenic neu- 
trinos can be translated into bounds on the fundamental Planck scale. ^3 
In the case of n extra spatial dimensions compactified on T" with a com- 
mon radius, the bounds derived using AGASA-|-Fly's Eye exposure represent 
the best existing limits on the scale of TeV-gravity for n > 4 extra spa- 
tial dimensions. A summary of the most stringent present bounds on AId 
for n > 2 extra dimensions is given in Fig. 13 Certainly, the lack of ob- 
served deeply-penetrating showers can be used to place more general, model- 
independent, bounds on cTiyAr.ES 



Up to now we have only discussed how to set bounds on physics beyond 
the SM. An actual discovery of new physics in cosmic rays is a tall order be- 
cause of large uncertainties associated with the depth of the first interaction in 
the atmosphere, and the experimental challenges of reconstructing cosmic air 
showers from partial information. However, a similar technique to that em- 
ployed in discriminating between photon and hadron showers can be applied 
to search for signatures of extra-dimensions. Specifically, if an anomalously 
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large quasi-horizontal deep shower rate is found, it may be ascribed to ei- 
ther an enhancement of the incoming neutrino flux, or an enhancement in 
the neutrino- nucleon cross section. However, these two possibilities may be 
distinguished by separately binning events which arrive at very small angles 
to the horizontal, the so-called "Earth-skimming" events.''^ An enhanced flux 
will increase both quasi-horizontal and Earth-skimming event rates, whereas 
a large BH cross section suppresses the latter, because the hadronic decay 
products of BH evaporation do not escape the Earth's crusti^ 

6 The PAO inquisition 

The Pierre Auger Observatory (PAO) is designed to work in hybrid mode, 
employing fluorescence detectors overlooking a ground array of deep water 
Cerenkov radiators. During clear, dark nights, events will be simultaneously 
observed by fluorescence light and particle detectors at ground level. The 
PAO is expected to measure the energy, arrival direction and primary species 
with unprecedented statistical precision. It will eventually consist of two sites, 
one in the Northern hemisphere and one in the Southern, each covering an 
area of 3000 km^ and each consisting of 1600 particle detectors overlooked 
by 4 fluorescence detectors. For showers with zenith angle < 60°, the overall 
acceptance (2 sites) is 14000 km^ sr. The angular and energy resolutions of 
the ground array (without coincident fluorescence data) are typically less than 
1.5° and less than 20%, respectively. "Golden events," those detected by both 
methods simultaneously, will have a directional reconstruction resolution of 
about 0.3° for energies near 10^^ GeV. If an event trigger is assumed to require 
5 detectors above threshold, the array is fully efficient at 10^" GeV. In three 
years of running, the surface arrays in both hemispheres will collect more than 
1000 showers above 10^°'^ GeV with approximately uniform sky exposure. 
This will enable a straightforward search for correlations with discrete sources 
and also a sensitive large scale anisotropy analysis 

For showers with zenith angle exceeding 60°, the aperture of PAO in- 
creases by roughly 50%. For a pure hadronic composition, there will be over 
1000 well reconstructed events beyond 10^" GeV, with a mean energy error 
~ 25%. On the other hand, if 7-rays are dominant at high energy, the rate 
will be reduced by an order of magnitude allowing for a clear discrimination 
between these two cases. PAO will be able to establish strong bounds on the 
7-ray flux at energies as high as 10^^ GeV.El 

In addition, PAO offers a window for neutrino astronomy above 10^ GeV. 
For standard neutrino interactions in the atmosphere, each site of PAO reaches 
^ 15 km'^ w.e. sr of target mass around 10^° GeV, which is comparable to other 
neutrino detectors being planned.'''^ Moreover, the sensitivity of PAO could 
be significantly enhanced by triggering on neutrinos that skim the Earth, 
traveling at low angles along chords with lengths of order their interaction 
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lengthM^As can be seen in Fig. [3 PAO will provide us with statistics to begin 
discriminating among the many promising ideas so far proposed to explain 
the origin and nature of CRs above the GZK energy limits. Measurements of 
quasi-horizontal neutrino fluxes will also allow better limits to be placed on 
low scale gravity (see Fig. 13 . 

7 Coda 

Statistics and better experimental handles should enable us to reconstruct 
the energy spectrum beyond 10^^ GeV, to locate the CR sources in the sky, 
and to discern the primary chemical composition. Future CR data will not 
only provide clues to the CR origin, but could enhance our understanding of 
fundamental particle physics. For example, if CR primaries are found to have 
a significant 7-ray component above 10^^ GeV, this could suggest an exotic 
ingredient in CRs, such as the decay products of TDs/MSRs, and thus could 
provide insight into the description of the early universe as well as particle 
physics beyond the SM. The puzzle of ultrahigh energy CRs may even have 
something to say about issues as fundamental as local Lorentz invariance: 
the absence of photo-pion production above the GZK-limit would imply no 
cosmogenic neutrino flux and possibly undeflected pointing of the primary 
back to its source. Additionally, contrasting the observed quasi-horizontal 
neutrino flux with the expected neutrino flux can help to constrain TeV-scale 
gravity interactions and improve current bounds on the fundamental Planck 
scale. An optimist might even imagine the discovery of microscopic BHs, the 
telltale signature of the universe's unseen dimensions. We are entering this 
new High Energy Physics era with the Pierre Auger Observatory. 
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